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Abstract. 

From the change in kinetic energy induced by an external field, we discuss 
the applicable conditions for the Mott-Davis and Moseley-Lukes form of the Kubo- 
Greenwood formula (KGF) for the electrical conductivity which has been implemented 
in ab initio codes. We show that the simplified KGF is suitable only for computing 
the AC conductivity at sufficiently high frequency and when the gradient of the carrier 
density is small. 
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1. Introduction 

The Kubo-Greenwood formula (KGF) is widely adopted to compute AC conductivity for 
amorphous semiconductors and to estimate the DC conductivity from an extrapolation 
procedure [H [2]. Greenwood[3] based his derivation on a kinetic expression for current 
density which is only justified when the gradient of the carrier concentration is smalljl]. 
Mott and DavispQ, Moseley and Lukes[5] developed a simplified version of KGF 
without using an explicit expression for the current density, but they made two implicit 
assumptions: (1) the Joule heat produced by a sample per unit time is f2j ■ E, where Q 
is the volume of sample, E the strength of external electric field, j is the current density; 
and (2) identifying the power Qj ■ E of current as the absorption energy per unit time 
r from AC field: 

F = Yl ^fi{ w fi p i ~ w if P f), (!) 
fi 

where Pi is the occupation probability of the initial state |z), Wfi is the transition 
probability per unit time from initial state \i) to final state |/), fkOfi is the energy 
difference between |/) and The conductivity is read off from 

7T a ^E a E p = T, a,p = x, y, z, (2) 

where the factor 1/2 in the L.H.S comes from averaging the power of current over one 
period of the AC field. The AC conductivity at frequency u is[U [5] 

= 2ne 2 tfn r dE [f(E) - f(E + hu)]\D\l N(E)N(E + jjoQ ^ 
m 2 J hjj 
where Q is the volume of sample, / is the Fermi distribution function, N(E) is the 
density of states, D E i >E = J d 3 xip E ,dipE/dx, and 'av' represents an average over all 
states having energy near E' = E + fko. The DC conductivity is obtained by taking the 
u -> limit in ©H]: 

**> = J dE\D\l v [N(E)] 2 ^- (4) 

Since Wfi contains an energy-conserving delta function 5(Ef — Ei — fiw), the assumption 
(2) means that the system is driven by a radiation field with frequency u = Ufi. The 
photon energy for a DC field is zero, so that the absorbed energy of the system from 
a DC field should vanish. According to ([2]), would be zero. However a DC voltage 
does produce Joule heat in a conductor or a semiconductor. One may conclude that 
(i) (EJ) is not suited to AC fields with low frequencies; and (ii) (j4j) is not consistent 
with (j2J), although the limiting procedure from (j3J) to (jlj) seems legitimate. On the 
other hand, ([3]) has been used to compute AC conductivity and extrapolate <Jd c in liquid 
carbon[2]. P| is implemented in SIESTA|6J to calculate adc for a-Si and a-Si:H. The 
DC conductivities obtained fall in the range of observed values in different systems. 
Therefore extrapolation from ([3]) and (jlj) must represent the correct adc to some extent. 
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The aim of this paper is to explore and resolve these controversial issues. In Secf2J 
we derive the rate dKj dt of change kinetic energy from the time-dependent Schrodinger 
equation. This rigorous expression of power is applicable to an arbitrary electromagnetic 
field, and leads to a proper current density [7J [8j 0]. Using the approximation implied by 
assumption (1), the rigorous current density reduces to the kinetic expression used by 
Greenwood[3]. Thus the simplified derivations [U E] suffer the same error as the original 
one: the formula is suitable only when the gradient of the carrier density is small, a 
condition derived from the current operator [I]. 

To find the connection between dK/dt (well defined for any electromagnetic field) 
and T which is only defined for a radiation field, we calculate dKj dt for three successively 
more complex cases: DC voltage, AC voltage and arbitrary field. For both DC and low 
frequency AC voltages (i.e. any field described by a scalar potential), dK/dt comes from 
the work done by the internal force. When a time-dependent vector potential appears, 
there is an additional dK/dt term contributing to the change in kinetic energy which 
results from the time dependence of vector potential dA/dt. 

A DC voltage described by a time-independent scalar potential can be described 
from both stationary perturbation theory (SPT) and time-dependent theory (TDPT). 
Secj3] will show that in TDPT only when we introduce the interaction with the 
time- independent scalar potential in a specified way (physically reasonable), are the 
descriptions of states and dK/dt in TDPT consistent with those in SPT. This is a 
natural requirement for a self-consistent theory. In SecJH we show that the Joule heat 
from a low frequency AC field described by a time-dependent scalar potential is the 
work done by the internal force. As expected, when u = 0, the result for AC voltage 
reduces to the DC result. This is possible because we consistently use a scalar potential 
to describe the interaction with a low frequency external field. If we adopt a different 
gauge, both scalar potential and vector potential are needed. In SecJH we calculate the 
rate of change of the kinetic energy induced by a general electromagnetic field described 
by both vector and scalar potentials. We will see that assumption (2) renders the 
simplified KGFpQ 5J ([3]) only suitable for computing conductivity at sufficiently high 
frequency, i.e for a radiation field. This is the origin of the confusion: viewing F as the 
unique reason for change in kinetic energy [U [5] dK/dt and taking the u> — > limit in 
(j3J)[lJ are not valid. Our formulation is applicable to an arbitrary electromagnetic field 
in any gauge: applying Greenwood's gauge A = — E£ for a DC voltage in the dK/dt 
term of our strict formula does lead to However (jlj) neglected the contribution from 
the work done by the internal force. 

2. Rate of change in kinetic energy 

We will use the Schrodinger picture. Consider a system with N electrons + Af nuclei in 
an external electromagnetic field (A, 0), at time t, and the state of system is described 
by ^'(r 1 ,r 2 , • ■ ■ ,r N ;t). To save space, we will not write out the nuclear coordinates 



4 



explicitly. \E' / satisfies the Schrodinger equation 

iffflf'/dt = H'W, (5) 

where H' = H + H m f(t), H m f(t) is the interaction between system and external field 
(A, 0). The time dependence of H m f(t) comes from the external field. H is the 
Hamiltonian of the system without external field: H\m) = e m \m), we use \m) or ty m 
to denote the m th stationary state of the system. If the system is in a thermal bath at 
temperature T, before introducing external field, the probability that the system is in 
state \l/ m is P m = e~ l3£m /Z, where Z = Y^ n e _/3en is the partition function. 

For a system in external field (A,0), the velocity operator v, for the i electron 
is[9] Vj = m -1 p™, where p™ = —iKV ri — eA(rj; t) and are the mechanical momentum 
and position operators of the i th electron, e = —1.6 x 10 _19 C is the charge of electron. 
Similarly v a = p™/M a is the velocity of the a th nucleus, p™ = -tW^ + Z a eA(R a ; t) 
and R Q are the mechanical momentum and position operators of the a th nucleus, —Z a e 
is the charge of the a th nucleus. The average kinetic energy of the system in state *&'(t) 
is 

Ky,(t) = J dT^'*(t)K(t)^'(t), (6) 
where dr = dridr', dr' = d?2 ■ ■ ■ dr^, the arguments of are (ri, rg, • • ■ , r^; t), and 

i a 

is the kinetic energy operator of the whole system in an external field. The time 
dependence in K(t) arises from that of A. To compute the macroscopic response to 
a mechanical perturbation, the coarse-grained average and ensemble average can be 
done in the final staged EJ E]. In this paper we only discuss the average over the state 
of the system. 

With the help of (JHJ), the rate of change in the average kinetic energy is [9] 

jKv{t) = 1 J dT*'*(t)[K,H')V(t) + J dr^*(t)^-^(t). (8) 

The Hamiltonian H' of system in external field can be written as 

H' = K{t) + V (j> + H 2 , (9) 

where 

V^ = ^e0(r,;t)-^Z Q e0(R a ;t), (10) 

i a 

is the potential energy of the system in an external field (A, </>), and 

ij ia a/3 

represents the internal interactions of the system, ^(r^Rj is the interaction energy 
of the i th electron and a th nucleus, V^r^r,,) is the interaction energy between the i th 
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electron and the j th electron, f/ 2 (R a , Rg) is the interaction energy between the a th and 
the (3 th nuclei. The I s * term in (JSJ) can be changed to: 

[K,H'] = [K,H 2 ] + [K,V li> ]. (12) 

Here (ih) -1 ^, H 2 ] represents the power due to the internal force: 

l[K , H 2 ] = • v, - » V„ • « + £«, • v„ - J-V^ • O, (13) 

i a 

where 

f Q = -{VrJ^V^R,) + ^[/ 2 (R a ,R^]}, (14) 

is the internal force on the a th nucleus, and 

S = -{V r J~ £ V 2 {v u vj) + £ Vi(r„ R a )]}, (15) 

j a 

is the internal force on the i th electron. From now on we will not write out the 
corresponding terms for nuclei which are similar to those for electrons. The second term 
of (fT2"|) represents the power due to the electric field described by the scalar potential: 

±_[K, Vy = e[-V ri 0(r i; t)} • v, + ^ ^^[V r , • V ri 0(r 4 ; t)}. (16) 

i i 

To calculate the 2 nd term in (JSJ), one should notice that [p m ,dp m /dt] ^ and 
[p- dp m /M}^0: 

Substituting ffTMBTT?]) into (JSJ, one finds 

^**'(*) = 1 dr*'*(*) ^ eE(r 4 ; i) • Vi *'(*) (18) 

i 

i i 

where 'E(r i ;t) = — V0(i"j;t) — dA(ri]t)/dt is the electric field at r 4 . To obtain the 2 nd 
and 4* h terms in (|18[) . we integrated by parts. (ITS]) is a form of the Ehrenfest theorem: 
the rate of change in kinetic energy equals the work done per unit time by the external 
electric field (the first two terms) and the internal force (the last two terms). The I s * 
and 3 rd terms are the corresponding quantum average values of the powers in classical 
mechanics. The 2 nd and 4 th terms are produced by the commutation relation between 
momentum and position: they will disappear in classical limit. 
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In the first two terms of (IT5|) . exchanging the integration variables rj,(k = 
2,3, ■ ■ • N) <H- r 1; using the antisymmetry of the many-electron wave function and 
changing ri to r, they become J n dr [E(r; £)] • j m (r;i): 

j m (r; t) = Nej dr'VvV + ^V r n'(r; t), (19) 

where the arguments of ^' are (r, r 2 r3 • • • rjy; t), n'(r;t) = N j dr'^f'^f'* is the carrier 
density j m defined by ( TT9T) is the same as the rigorous microscopic current density 
obtained from the microscopic response method[7] and polarization density [8J. Only 
when the gradient of the carrier density is small, can one neglect the second term 
in (fT9~j) and replace j n dr E • j m with Oj • E. By means of the equivalence between the 
microscopic response method and the Kubo formula [1], the first term in (fT9j) corresponds 
to the kinetic expression eTr[p'(i)v] of Greenwood[3], where p'(t) is the density matrix 
of the system in external field. Thus assumption (1) is equivalent to using eTr[p'(i)v], 
the kinetic expression for the current density. 



3. DC voltage 

Let us adopt the common gauge in which a DC field is solely described by a time- 
independent scalar potential. In this gauge, both stationary perturbation theory (SPT) 
and time- dependent perturbation theory (TDPT) can be used, and they should give the 
same results for any observable quantities. As we will see, the first nonzero contribution 
to the rate of change of kinetic energy is second order in V (consistent with our 
macroscopic experience), so to formulate a consistent approximation, we will carry out 
perturbation theory to second order in V. 

Applying a DC voltage on a piece of conductor or semiconductor, after a short 
transient period, the system will evolve to a steady state if the system is in good 
thermal contact with the environment such that the Joule heat evolved can be completely 
removed from the system. The constant external voltage establishes a time-independent 
electric field inside the system. The system is described by a Hamiltonian H' = H + V, 
where H is the Hamiltonian of system without external DC voltage, V is given by (ITUI) 
with a steady scalar potential 0. 

If the system is initially in an eigenstate tyj of H with eigenvalue Ej, after the short 
transient period, the system will be in the stationary state of H' with eigenvalue e'j. 
Denote Vkj = (^k\V\^/j), then one can easily compute[9] and e'j to second order in 
V: 

vI/ ; = ^+E{E c S VI/ -+ c S ) ^}- (2°) 

P=l,2 m&j) 

Since H' is time-independent, the time evolution of system is given by 

y'.(t) = < ' 7 ' h %. (21) 



Substituting \E^- and e'j obtained from SPT into (12"T1) . to second order in V, 

= e-^i/A^. + e"^ /ri ^ { J2 6 S*m + (22) 

p=l,2 m(=£j) 



where 



and 



and 



6 (2) _ KnfcVfcj f 2 4l 

T T V m j Zt V m j . 

(Ej - E m ) n Ej - E m 



,(2) _ _ L \- \V kj \ z it \- I Vfcjl^ _ f2 Vj§ / 9 ,n 

fc(^i) 1 J fej J 
Now consider the viewpoint of TDPT, in which the perturbation V causes 
transitions from ^ to other eigenstates of if. Using the familiar expansion[9j 

^.(t) = ^. e -^/A + ^ [ ^ a^(t)e- ite -/^ m + ag ) (t)e-^'/ fi ^], (26) 

p=l,2 m(^j) 

the first order expansion coefficients satisfy 

ih m J t = e it{£m -^ )/R V mj for m ^ j, (27) 

and 

da$(t) 

ih " = tt,-. (28) 
dt 13 

To make (1261) consistent with (122]) at order V, we have to integrate (1271) by adiabatically 
introducing the interaction dt'e xt ' (A — > + ) and integrate (1281 by dt'. This 
procedure is reasonable because the system is initially in state tyj, the transition from 
tyj to another state *S? m (m ^ j) requires some time. On the other hand, the probability 
amplitude of state tyj begins to decrease immediately. 
The second order expansion coefficients satisfy 

dn^ (t\ 

= £ a^{t)^-^V mk + <$(t)<*«"-<M*V mj , (29) 



and 
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If we integrate ( 121?j) by adiabatically introducing interaction f_ dt?e xt ' (A — > + ), 
integrate ( 130]) by J Q * <it' and in a?- drop one term with wrong time factor e~ tt£k ^ h } we 
almost reproduce f l22l) except without the 1/2 factor in the first term of . 

There are two differences between SPT and TDPT of concern to us here. In SPT 
no equation exists to determine cf^{p = 1, 2) in ( 1201) . The perturbed wave function (1201) 
is not normalized, if one does not include Ylp=i2 c jj^j- c jj(P =1)2) are determined 
from the normalization of to the corresponding order [9]. In TDPT, df){t) (p — 1, 2) 
are determined by (1281130]) . Using ©, one can easily find (d/dt)[f dr^ = 0. 
Thus the perturbed wave function ( 126]) is normalized if J dr^^j = 1. If in ( 126]) . we 
only considered o^j(t) for m ^ j, and used the normalization condition of ty'At) to 
determine afj(t) (p = 1,2), we would not reproduce ( 122]) . Therefore the suggested 
means of introducing the interaction is necessary to make TDPT consistent with SPT. 

When the DC field is described by a time-independent scalar potential, the rate of 
change in kinetic energy can be written as 

jK %{ t) = drV;{t)[K^H]%{t) drV*{t)[H,V]%{t), (31) 

where K is the kinetic energy when vector potential is zero. Because [K , H] = [K , H 2 ] 
and [H, V] = [K , V], after a comparison with (I13|ll6p . one may say that the first term 
in ( 13~TT) is the power due to the internal force, and the second term in ( 13TT) is the power 
due to the external force. Of course, the effect of external field is also reflected in ^-(t)- 
With the help of ff26|) . one can easily show that to order V 2 , the second term in f l3~T|) is 
zero. ( 13T1) becomes 

jK^t) = ^f drV*(t){K , H]%{t) (32) 
= k a kj* a< ifl £ i - £ k]K ki + rr J2i £ i ~ £ j][ a( if K oji ~ a \f* K oijl 

kl I 

where K 0j i = (j\K \l) and off = af^ (t)e^ ej ~ ei ' }t / h (p = 1,2). The change in kinetic 
energy is produced by the power of the internal force. It is easy to check that 
without an external field, the internal force does no work, the average kinetic energy 
K^ d (t) = J dT^/*(t)Ko^j(t) in stationary state *&j(t) does not change with time: 

jK^{t) = ±_J dT**(t)[K ,H]%(t) = 0. (33) 
The Joule heat comes from a steady voltage which changes the state of system. 

4. Low frequency AC voltage 

We will use a gauge in which a low frequency AC field is solely described by a time- 
dependent scalar potential. The interaction of system with an AC voltage is given by 
ffTUj) with a periodic scalar potential 0: 

V(t) = Fe~ iwt + (34) 
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An AC voltage will produce a time-dependent current. According to Ampere's law, the 
time-dependent current will produce a time-dependent magnetic induction. Therefore 
a time-dependent vector potential A(r, t) must accompany the AC voltage. ( |34|) is only 
suitable for low frequency u <C cr/e , where a is the DC conductivity for the system. 

When the frequency of an AC voltage approaches zero, its properties should be the 
same as those of DC voltage. Therefore we must integrate the equations of probability 
amplitudes cbfjit) (p = 1, 2) for V(t) in the same way as those for the DC voltage: 



^ l '(f) = lj ill 

'•' ' h(uij -u) h(uij + u) 



«w (*) = -tF — - \ for I ? j, (35) 

and 



FjA e-** - 1) _ - 1) 



^M = W .."' -^^. (36) 



We will not write down the expressions for a\- it) and a^- (t), they are too long. To the 
2 nd order in V(t), the state of system is 

ttj(t) = V je - it£ i /h a\f(t)^ie~ it£l/h . (37) 

p=l,2 I 

The rate of change in kinetic energy is 

j f K % (t) = dr%*(t)[K ,H}%(t)+^ J dr%*(t)[H,V(t)}%(t),(S8) 

For an AC voltage, to obtain the dissipated energy, we must average (138]) over a period 
T = 2-k/u of the AC voltage [10]: T" 1 jf T dt. Using a\f{t) and af]{t) (p = 1, 2), one can 
show that to 2 nd order in V(t), 

T- 1 [ dt\ [ dT*'j , {t)[H,V(t)]*' j {t) = 0. (39) 
Jo i>n> J 

The time-averaged power of the external force is zero. Thus the change in kinetic energy 
is due to the internal force: 

t So dt Jt K ^ {t) = rl Tdt k J ^;w[#o,#wo (40) 

_ J_ v-^ r K \- r F mk F* k F km F kj 1 (FmjFjj - F* m Fjj)K jm 



_^ r Frnk^jk F km F kj ^ (F^jFjj — FjmFj^Kmj 

1713 / -— ' ^CJfcj + Cj) ^(Wfcj — to) huj 



fcj(^n h(uj kj - Ujh^UJij - U) h{u kj + UJjh^UJij + LO) 

For an AC voltage, because [K Q ,H] = [K Q ,H 2 ], we may say that the Joule heat comes 
from the power of the internal force. For u> = 0, Eq. (jlO"l) reduces to Eq. (l32l . the Joule 
heat for DC voltage. 
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5. Electromagnetic field 

In Secj3]and|H a special gauge is used: both DC voltage and low frequency AC voltage 
are described by scalar potentials. The rates of change kinetic energy are given in 
( !3TTl38l) . Now consider the system interacting with a general electromagnetic field 
described by (A, 0) which may or may not change with time. We will not restrict 
ourselves to any special gauge. The kinetic energy operator of the system is K(t) 
rather than Kq. If vector potential A changes with time, there is one more term 
dK{t)/dt in dK^.(t)/dt (jSJ). dK{t)/dt results from the time dependence of vector 
potential <9A(r, t) /dt. 

To apply TDPT to compute dK^(t)/dt, we notice that H m f = Va + V^, where 
Va = Vai + Va2 is the interaction involving vector potential A, Vai presents the terms 
which are first order in A, Va2 presents the terms which are second order in A. Vai 
is only a function of coordinates and does not include differential operators. Notice 
K = Kq + Va and H = Kq + H2, the commutator in the first term of (jSJ) can be 
transformed to 

[K, H') = [K, H] + [H, H mf ] + [Vai, V], (41) 
where V = V^ + H 2 and 

^Vai,V] = Y,rn~ 1 [eA(r t ,t)-V r V]-J2M- 1 [Z a eA(R a ,t)-Vn a V].(42) 

i a 

Since — V r V = f» + [— eV ri 0(r i; t)], — eA(r^ t)/m is the part of velocity due to field of 
the i th electron, the first term in ( )42i) is the power of the field momentum due to scalar 
potential and the internal force. The 2 nd term is the power of the field momentum of 
nuclei. Substituting ( )4T1) into (|E|), one has 

jK*>{t) = iy dT%*(t)([K,H] + [H,H mf (t)])^(t) (43) 

+ jdr^(t)(^ + l[V A1 ,VW ] (t). 

For an electromagnetic field with several frequencies u n , the matrix element of the 
interaction has form: 

[Hmf{t)]l 3 = Y} F Z e ~ Wnt + iT e ***]- ( 44 ) 

n 

With the same method for an AC voltage described by a scalar potential, one can show 
to 2 nd order in H mf (t), (ih)- 1 J dr^'*{t)[H, H mf )^'(t) = 0. The first term in flED can 
be similarly obtained as fj40l for AC voltage: 

dTV;{t)[K,H\%{t) = ^ E^Efe-^f We^^' )/ %-(45) 

P=1.2 
kl{k^l) 
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where the matrix elements of the kinetic energy are calculated with K(t) rather than 
Kq, the transition probability amplitudes are computed for H m f(t). After average over 
one period of external field, the order V term in the curly bracket is zero. We see from 
(jT7|) that dK{t)/dt is first order in vector potential. To second order in H m f(t), 

/** ? ( t) ^*; W = /«-^ (46 ) 

+ 5>£ ) *(i)e ,,(! <-'> / « / + ag ) (i)e-" < *'-* j)/S J drt'^-n- 

For many choices of gauge, dK{t)/dt is Hermitian, the second term in the square bracket 
is the complex conjugate of the first term. Combining (j42|45|46p . the rate dKy.(t)/dt 
of change in kinetic energy in ( |4~3]) is determined. 

We analyze the conditions which lead to KGF for this general case. If (i) the 
gradient of the carrier density is small; and (ii) the wavelength of vector potential is 
larger than the characteristic length of the considered sample, the first term in (jJEJ) is 
zero. One can see this from ffl7|) : under conditions (i) and (ii), the second term in ffTTl) 
is ignored, and the first term in ( l4U|) becomes [— 9A(r; t)/dt] J2i I dr^jV^j. But the 
average velocity J G^r\F*Vj\F,- in a stationary state tyj of H is zero. Now only the square 
bracket term is left in (H6|) . Each term represents the power due to part of electric field 
described by vector potential: 

- dA(r;t)/dt = -dA ± (r;t)/dt - dA\\(r;t)/dt, (47) 

the longitudinal and transverse parts satisfy V x An = and V • Aj_ = respectively. 
The contribution from — V0 is absent from (1461) . 

Greenwood used a special gauge: A(t) = — tE and (p = to describe a DC 
voltage[3]. Then V A1 = tE ■ v, according to TDPT, (l\eEt ■ v|j) ~ (e, - e^af) /t. 
If the interaction time is long enough that transition probability per unit time is well 
defined, then a typical term in the square bracket of fT46l) becomes J2i( £ i ~ £ j)\ a tj > \ 2 1^- 
Averaging (1461) over the occupation probability of the initial state tyj, one obtains (jl]). 
From (I18|H9D and (T43]), (T4]) missed the contributions from (S2J) and (J45]). Although the 
momentum due to field i.e. ( 1421) is negligible, the contribution f j4"5l) is same order as 

m- \(VAih\ 2 /h. 

If we consider only the radiation field in (l47j) . dA±(r;t)/dt ~ uA±. Then 
a[f(t) = -(/|eA_L • v|j)e i( ^ _w) ' +A *[/i(w; j - u - zA)]" 1 . If the interaction time is long 
enough, the square bracket in (|4"6"1) becomes hw\a^\t)\ 2 /t i.e. ([T]). For a zero 
frequency radiation field, the absorbed energy from field is zero. Thus we understand 
that although the limit procedure from Q to (jlj) is not legitimate, (jl]) can be obtained 
in Greenwood gauge (a longitudinal field), and contains part of the conductivity. The 
power fj45l) induced by the internal force on system always exists, no matter longitudinal 
field or transverse field. The applicable lower limit frequency of is at least to » cr/e . 
For intrinsic Si[TT], a = 1.2 x 10 _5 Scm _1 , the simplified version[TJ[5] of Kubo-Greenwood 
formula wok only when the frequency of external field is higher than cx/eo ~ 1.4 x 10 8 Hz. 
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Numerical caculations[2J H2J EH EEH EE] have shown cr(u;) approaches zero with 
uj below a resonance absorption frequency u r which is roughly the splitting from the 
maximum of density of states in the valence band to the edge of conduction band. This 
agrees with our analysis. Sometimes researchers extropolate[2] from a(cj r ). Because 
of the existence of (jlj), if the velocity matrix elements are not very sensitive to the 
energies of states, extrapolation will deliver reasonable results. 

6. Conclusion 

In summary, from the rate of change in kinetic energy, we obtained a rigorous expression 
for the power done by an arbitrary electromagnetic field in any gauge. It leads to a 
proper current density which has been proved by the continuity equation]?], polarization 
density [8] and current operator [4]. We show the simplified derivation of KGF by Mott- 
Davis and Moseley-Lukes suffers the same approximations used by Greenwood. The 
work done by the internal force is missed in (131141) . they are same order as the terms 
in KGF. Using (jlj) or extrapolation from (J3J) can obtain a significant part of the DC 
conductivity, but a stricter calculation based on the rigorous current density would 
deliver more accurate, possibly even qualitatively new results. 
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